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Abstract 

Using the quark-diquark approximation in the framework of Buchmiiller-Tye potential 
model, we investigate the spectroscopy of doubly charmed baryons: and . Our results 
include the masses, parameters of radial wave functions of states with the different excita- 
tions of both diquark and light quark-diquark system. We calculate the values of fine and 
hyperfine splittings of these levels and discuss some new features, connected to the identity 
of heavy quarks, in the dynamics of hadronic and radiative transitions between the states of 
these baryons. 
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1 Introduction 



For experiments at the running and planned high-energy machines such as LHC, B-factories 
(HERA-B, for example) and Tevatron with high luminosities, there are some hopes and pos- 
sibilities to observe and study new hadrons, containing two heavy quarks, like the doubly 
charmed baryons (ccq) (here and throughout this paper, q denotes a light quark u or d). 

Looking for such projects, it is important to have reliable theoretical predictions as a guide 
to the experimental searches for these baryons. In our previous papers we have already studied 
both the production mechanisms of doubly charmed baryons on the different future facilities 
fl]] , [|| and the total inclusive lifetimes of such systems || . In this work we would like to address 
the question on the spectroscopy for these baryons in the framework of potential approach in 
the Buchmiiller-Tye model. The (QQ'q) spectroscopy was also considered in some other models 



The baryons with two heavy quarks include the features of both the dynamics in the D- 
meson with a fast moving light quark, surrounding a static 3-color core, and the dynamics of 
heavy-heavy mesons (J/^, B c , T), with two heavy quarks sensitive to the QCD potential at 
short distances. So, a rich spectrum is expected. First, there are excitations due to the relative 
motion of two heavy quarks in the quark-diquark approximation, which we use here. Some 
comments on the validity of this approximation will be given later. Second, we consider also 
the excitations of the light quark or combined excitations of both degrees of freedom. 

As we will show in this paper, some new features, related with the identity of two heavy 
quarks, arise in the dynamics of doubly charmed baryons (or doubly beauty baryons). The 
hadronic or electromagnetic transition between the 2P — > IS diquark levels is forbidden in 
the absence of interaction of diquark with the light quark or without taking into account 
nonperturbative effects. The qualitative picture of this effect will be discussed below, and a 
rigorous quantitative solution of this problem will be given in our subsequent papers. 

This work is organized as follows. Section 2 is devoted to the determination of masses, 
values of radial wave functions (for S'-levels) and their first derivatives (for P-levels) at the 
origin for both the diquark and the system of light quark-diquark. In Section 3 we evaluate 
the fine and hyperfine splitting of different levels for the above systems. Section 4 contains 
some comments on the radiative and hadronic decays of doubly charmed baryons. And, finally, 
Section 5 draws our conclusion. 



Investigating the baryon spectroscopy, one faces a three body problem in the framework of 
nonrelativistic quantum mechanics. Its reduced hamiltonian has the following expression: 



2 Mass spectrum of doubly charmed baryons. 





where x, y are Jacobi variables: 



x = r 2 -ri, 




y 




(3) 



where M is the heavy quark mass, and m is the light quark mass. 
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There are several methods for a solution of three-body Schrodinger equation ||. The 
first way is the variational methods, where one expands the wave-function in terms of the 
eigenstates for a symmetric harmonic oscillator, or gaussians, at different ranges in the x and 
y coordinates, or uses the hyperspherical formalism. The other methods are the quark-diquark 
approximation and Born-Oppenheimer approximation. The former is used in our evaluation 
of doubly charmed baryon spectrum. The reason is that the ground state of (ccq) consists of 
a localized (cc) cluster surrounded by the light quark q, with the average distance (r cc ) much 
less than {r cq ). However, when one accounts for the radial or orbital excitations of diquark, 
the average separation between the heavy quarks increases, and the quark-diquark structure 
is distroied. So, in this region, our results for the mass spectrum of these baryons are quite 
crude. 

Next, a dramatic simplification of (ccq) dynamics is obtained in the Born-Oppenheimer or 
adiabatic approximation. Two heavy quarks have much less velocity than that of the light 
quark. When they move, the light quark wave function readjusts itself almost immediately 
to the state of minimal energy. Therefore, the calculation can be done in two steps: for any 
given x, one computes the binding energy e(x), which is, then, used as an effective potential 
governing the relative motion of heavy quarks. This was done in GJ: first, in the nonrelativistic 
potential model and, second, in a variant of MIT bag model. From our point of view, this 
method is the most suitable for the baryon spectroscopy. However, in this work, we will be 
satisfied with the accuracy, given by the quark-diquark approximation. 

In present work, we use the QCD-motivated potential || given by Buchmuller and Tye ||, 
which was justified on the J j^f and T spectra with the following values of parameters: 

m c = 1.486 GeV, m q = 0.385 GeV, (4) 

where the light quark mass is obtained by the fitting of theoretical predictions with the 
Buchmuller- Tye potential for the Z)-meson mass to its experimental value. Of course the 
motion of light quark is relativistic inside the D meson, as well as inside the baryons under 
consideraton and so can not be treated in the framework of nonrelativistic quantum mecanics. 
But we think, that considering its motion in both cases in the same way, we can obtain good 
approximation for the values of mass levels of our baryons. 

Estimating the mass spectrum of diquark excitations, one has to take into account the 
factor of 1/2 for the potential due to the fact, that the heavy quarks inside the diquark are in 
color antitriplet state. At distances ~ 0.6 — 0.8 fm we have an attraction of quarks inside the 
diquark, so we assume that the shape of potential in this region related with their pairwise 
interaction is the same as for quark - antiquark interaction in heavy mesons. At larger distances 
we can not tell anything about it, so our estimates for higher-lying levels is quite rough. But 
for low-lying energy levels in this system the wave function is already zero in the region of 
large distances, so that our predictions for them can be trusted. Then, solving the Schrodinger 
equations for the diquark and light quark excitations, one finds the diquark and E+- 
baryons mass spectra and the characteristics of radial wave functions for both the diquark and 
the system of light quark-diquark Rd(0), Ri(0), R'd(9) an d R'i(0), shown in Tables 1, 2 and 3. 

In calculations we assume, that the threshold mass value of doubly charmed baryons is 
determined by the hadronic decay into A c -baryon and .D-meson, and, hence, it equals 4.26 



GeV [10]. The threshold for the stability of diquark is estimated from the following result, 



stated for a heavy quark-antiquark pair [11|: if a heavy quark and a corresponding antiquark 
are separated by a distance greater than 1.4-1.5 fm, then the most energetically favorable 
and probable configuration results in a production of light quark-antiquark pair, which leads 
to the fragmentation into a pair of flavored heavy mesons. So, we suppose that the same critial 
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Diquark state 


Mass (GeV) 


(r 2 ) 1 / 2 (fm) 


Diquark state 


Mass (Gev) 


(r 2 ) 1 / 2 (fm) 


IS 


3.16 


0.58 


3P 


3.66 


1.36 


2S 


3.50 


1.12 


4P 


3.90 


1.86 


3S 


3.76 


1.58 


3D 


3.56 


1.13 


2P 


3.39 


0.88 


4D 


3.80 


1.59 



Table 1: The (cc)-diquark spectrum: masses and mean-square radii. 

distance scale can be used for the colored diquark system, which results in the fragmentation 
of diquark to the heavy meson and heavy-light diquark. 



3 Spin-dependent splitting. 

In accordance with the results of refs. [h|] and [[14] , we introduce the additional term to the 
potential to take into account the spin-orbital and spin-spin interactions, causing the splitting 
of nL-levels in both the diquark and light quark-diquark system (n is the principal quantum 
number, L is the orbital momentum). So, it has the form 



2 V2m 2 



dV{r) 8 1 



rdr 



3 
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for the diquark splitting, and 
1 /L • S d 
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1 



x [6(L-S) 2 + 3(L-S) -2L 2 S 2 



4L 2 



-6(L-S d ) 2 -3(L-S d ) + 2L 2 S d 2 ])^ 



1 



n (diquark) 


R din s)(0)(R' d(nP) (0)) 


n (diquark) 


^(nS)(0)(^ fnP) (0)) 


IS 


0.530 


2S 


-0.452 


2P 


0.128 


3P 


-0.158 



Table 2: The characteristics of diquark radial wave functions Rd( n s)(Q) (in GeV 3 / 2 ), R'd( n p)(fy ( m 
GeF 5 / 2 ). 
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nd(diquark) 
ni (light quark) 


Mass (GeV) 


(*W°)) 


rid(diquark) 
ni(light quark) 


Mass (GeV) 


R l{nS) (0) 

(^W°)) 


IS IS 


3.56 


0.499 


IS 2P 


4.03 


0.118 


2S IS 


3.90 


0.502 


2S 2P 


4.36 


0.119 


3S IS 


4.16 


0.505 


3S 2P 


4.62 


0.121 


2P IS 


3.79 


0.501 


2P 2P 


4.25 


0.119 


3P IS 


4.06 


0.504 


3P 2P 


4.52 


0.119 


3D IS 


3.96 


0.503 


3D 2P 


4.42 


0.117 



Table 3: The mass spectra and characteristics of light quark radial wave functions in the doubly 
charmed baryons: S+ + and S+ with the different excitations of diquark, rid, an d light quark-diquark 
system, n/. 

for the light quark-diquark system, where V(r) is the phenomenological potential (Buchmiiller- 
Tye (BT) potential in our case), Si and Sd are the light quark and diquark spins, respectively. 
The first term in both expressions takes into account the relativistic corrections to the potential 
V(r). The second, third and fourth terms are the relativistic corrections coming from the 
account for the one gluon exchange between the quarks. a s is the effective constant of quark- 
gluon interactions inside the baryons under consideration. 

Expression (6) for the additional part of potential, causing the splitting of levels in the 
light quark-diquark system is obtained from the summing of the pair interactions, like (5), for 
the light quark with each of the heavy quarks. We include also the correction, connected to 
the interaction of internal diquark orbital momentum with the light quark spin. 

The value of a s parameter in (5), (6) can be determined in the following way. The splitting 
of the S"-wave heavy quarkonium (Q1Q2) is determined by the expression 

AM(nS) = la s ^—\R nS (0)\ 2 , (7) 
9 m\mi 

where R n s(fy is the radial wave function of the quarkonium, at the origin. Using the experi- 



mental value of liS-state splitting in the cc system [15| 



AM (IS, cc) = 117 ± 2MeV (8) 

and the Ris(0) value calculated in the potential model with the BT potential for the cc system, 
one gets the value of a s (^>) coupling constant for the effective Coulomb interaction of heavy 
quarks. 

In the present paper, we take into account the variation of the effective Coulomb coupling 
constant versus the reduced mass of the system (jj). In the one- loop approximation at the 
momentum scale p 2 , the 'running' coupling constant in QCD is determined by the expression 

where 6 = 11 — 2nj/3, and rif = 3, when one takes into account the contribution by the virtual 
light quarks, p 2 < m 2 . We assume that the average kinetic energies of c-quarks inside the 
diquark and the light quark-diquark system (the average kinetic energy weakly depends on the 
reduced mass of the system) are equal to: 

(T d ) ~ 0.2 GeV, (10) 
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(71) « 0.4 GeV, (11) 
so that, using the expression for the kinetic energy, 

(T) = ^ (12) 
where fi is the reduced mass of the system, one gets 

a ' (p2) = fc.ln(2<rWA &CD )' <13) 

so that a s (#) = 0.44 at Aqcd ~ 113 MeF. 

As one can see from equations (5) and (6), in contrast to the LS'-coupling in the diquark, 
there is the jj-coupling in the light quark-diquark system, where diquark and light quark have 
the different masses (here LS/ is diagonal at the given J; momentum, (J; = L + S/, J = J/ + J), 
J is the total spin of the system, and J is the total spin of the diquark (as we will see below 
in the case of interest, J equals to S^)). 

To calculate the values of level shifts, appearing because of the spin-spin and spin-orbital 
interactions, one has to average expressions (5), (6) over the wave functions of the correspond- 
ing states. Then, because the leading contribution to the spin-orbital splitting of the light 
quark-diquark system is given by the term + fcKsfj), we can use the state vectors 

with the given values of J and J; as the first approximation to the eigenvectors of the potential. 
For the potential terms, which are not diagonal in these states, we can choose another basis of 
vectors with the given values of J and S = S; + J 



| J; Ji) = Y,(-l) (J+Sl+L+J) ^/(2S + 1)(2J, + 1) \ '[ b j * \ \J;S) (14) 

or J and Jd 




\J-,Ji) = ^(-l) (J+Sl+L+J W(2J d + l)(2J l + l)^ * ^ j^l^Jd) (15) 

so that the potential terms of the order of l/m c m/, 1/m^ lead, generally speaking, to the 
mixing of levels with the different Ji values at the given J values. The identity of heavy quarks 
results in Sd = 1 at even and = at odd Lj, where is the diquark orbital momentum. 
Table 3 shows us that we has to take into account only the spin-orbital splitting of 1S2P and 
3D IS levels. 

In the first case (the splitting of the light quark-diquark system levels A^ J ) for 1S2P) one 
has: 

A ( l } = 17.4 MeV. (16) 

The levels with J = | (or i), but with the different J;, get the mixing. For J = |, the mixing 
matrix is equal to 



4.3 -1.7 
-1.7 7.8 



MeV (17) 



with the eigenvectors 



|152P(|')) = 0.9861^ = ^ + 0.1641^ = 1), (18) 
\IS2P(^)) = -0.164|J, = |> +0.986| J, = i>, 
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and the eigenvalues 



A; = 3.6 MeV, (19) 
Ai = 8.5 MeV. 



For J = \, the mixing matrix equals 

( -3.6 -55.0 



with the eigenvectors 



and the eigenvalues 



-55.0 -73.0 / MeV ™ 



|1S2P(1')> = 0.957| J", = |> - 0.291| Ji = |), (21) 
|152P(i)) = 0.291| J", = |> + 0.957| Ji = 



A' 2 = 26.8 MeV, (22) 
A 2 = -103.3 MeV 



In the second case (the splitting of diquark levels A( Jd ) for 3-D1S*) one gets: 

A( 3 ) = -3.02 MeV, 

A (2) = 2.19 MeV, (23) 
A (1) = 3.39 MeV 

(24) 

For the spin-spin interactions inside the diquark, we have only a shift of the energy levels 
because of the identity of heavy quarks. 

The hyperfine splitting for the light quark-diquark system can be computed, making use 
of the following formula: 



A«. = l(S(S + 1) - J(J + 1) - \)a s {2nT)^— * |i? z (0)| 2 , (25) 
' < 4 rn c mi 



where Ri(0) is the value of radial wave function in the system of light quark-diquark at the 
origin, and 

A?, = la s (2f,T)^*\R d (0)\ 2 , (26) 



for the diquark level shifts, where Rd(0) is the radial wave function of diquark at the origin. 

For the IS" and 2S"-wave states of diquark we have the following values of the shifts (because 
of the identity of the quarks there are no splittings of these levels) 

A(15) = 6.3 MeV, 
A(25) = 4.6 MeV 

The mass spectrum of doubly charmed baryons (H++ and £+), with the account for the 
calculated splittings is shown in Fig.l and Table 4. 
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Figure 1: The spectrum of doubly charmed baryons: H+ + and S+. 



4 Transition between diquark levels: 2P 
laboratory for the long distance QCD? 



IS 



Our plans for a future include a detail investigation of hadronic and radiative transitions in 
the spectrum of doubly charmed baryons. However, here we would like to discuss some new 
futures in the dynamics of baryons, containing two identical heavy quarks. 

So, because of the identity of two heavy quarks, we have a metastable 2P-wave diquark 
state. This state has the L = 1, S = quantum numbers, and, hence, the transition to 
the ground state (L = 0, S = 1) would require the simultaneous change of orbital and spin 
quantum numbers. It is worth to stress that the existence of such state is possible only in the 
baryons with two heavy quarks, because in ordinary heavy baryons the light diquark is never 
excited because of extremely large size. We have two possible scenario for the realization of 
such transition: 

1. A three-particle interaction via three-gluon vertex. This interaction breaks down quark- 
diquark picture and leads to some new wave functions of these states in the form: C\\L = 
1,S = 0} + C2I-L = 0, S = 1), where | C2 1 -C |Ci| for 2P. So, as one can easily see, the account 
for such interactions makes the radiative Ml-transition to the ground state to be possible. 
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(rid(diquark) 


Mass (GeV) 


(nd(diquark) 


Mass (GeV) 


ni(light quark)), J p 




niilight quark)), J p 




(IS lS)l/2+ 


3.478 


(3P lS)l/2~ 


3.972 


(IS lS)3/2+ 


3.61 


(3D lS)3/2'+ 


4.007 


(2P lS)l/2" 


3.702 


(IS 2P)3/2'- 


4.034 


(3D lS)5/2+ 


3.781 


(1S 2P)3/2~ 


4.039 


(2S lS)l/2+ 


3.812 


(IS 2P)5/2- 


4.047 


(3D lS)3/2+ 


3.83 


(3D lS)5/2'+ 


4.05 


(2P lS)3/2" 


3.834 


(IS 2P)l/2'" 


4.052 


(3D lS)l/2+ 


3.875 


(3S lS)l/2+ 


4.072 


(IS 2P)l/2~ 


3.927 


(3D lS)7/2+ 


4.089 


(2S lS)3/2+ 


3.944 


(3P lS)3/2~ 


4.104 



Table 4: The mass spectra and characteristics of light quark radial wave functions of doubly charmed 
baryons: S+ + and 5+ with the different excitations of diquark, n d , and light quark-diquark system, 

m. 

2. A nonpcrturbative transition given by an operator, proportional to [is r ■ V H(0) ■ (S\ — 
S2), where (jlb is the Bohr magneton, H(r) is the chromomagnetic field and r is the distance 
between two heavy quarks (here we have the same situation as in the case of transition between 
the orto- and para- hydrogen) . This transition goes with the emission of pion. We would like to 
stress the nonperturbative nature of this transition, because for its realization, the necessary 
condition is that the chromomagnetic field at the different points is correlated. The absence of 
such correlation prevents this scenario of transition from realization, because it would require 
two consequent gluon exchanges with the light quark or quark-gluon sea, wherein the orbital 
or spin quantum numbers of state would change, what is impossible because of the identity of 
two heavy quarks. 

A detail quantitative investigation of the nature for this transition is a subject of our sub- 
sequent papers. It will allow us to gain more information on the nonperturbative dynamics 
of QCD (in the case of second type transition) and, particularly, on the behaviour of chro- 
momagnetic field at large distances. So, we close this section with the question, posed in its 
title. 



5 Conclusion 

In this work we have used the Buchmiiller-Tye potential model within the quark-diquark 
approximation to describe the mass spectrum of doubly charmed baryons: H+ + and £+, 
including the fine and hyperfine splittings of mass levels. We have discussed the uncertainties 
involved and some possible ways to reduce them: the use of Born-Oppenheimer approximation 
and the corrected potential for the level splitting. In the previous section of article we have 
considered a new phenomena, arising in the radiative or hadronic transitions of 2P — > IS 
diquark levels. We have commented on some possible scenaria of this transition and lessons, 
we could learn from a studying the nature of such transitions. 
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